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Abstract 

We calculate inclusive photoproduction cross-sections for both proton and nuclei targets, 
using a saturation model derived from an approximate solution to the Balitsky-Kovchegov 
nonlinear evolution equation. This paper extends our hypothesis, that one can successfully 
use the dipole picture and saturation models to describe soft processes, instead of the soft 
Pomeron Regge parametrization. Our fit is compatible with our previous fit to DIS data, 
and utilizes the same phenomenological parameters as in our paper devoted to soft hadronic 
interactions. Using the Glauber formalism, and no additional parameters, we calculate 
the cross-sections for photoproduction on various nuclei, and compare our results with the 
relevant data. 



* e-mail 
t e-mail 

* e-mail 
§ e-mail 
^e-mail 



:jochen. bartels@desy.de 
gotsman@post.tau. ac.il 
leving@post.tau.ac.il; levin@mail.desy.de 
lublinm@mail.desy.de; mal@tx. technion.ac.il 

maor@post.tau.ac.il 



1 Introduction 



Over the past few years there has been much progress in the successful application of perturbative 
QCD to DIS processes for values of Q 2 > 2 GeV 2 . For very small Q 2 , Regge theory (e.g.[l]) 
provides a reasonable description of the data. Attempts have been made to bridge the two Q 2 
regions, using methods incorporating high parton densities and saturation. Models based upon 
these ideas have lead to excellent descriptions of the DIS cross-section for all values of Q 2 and 
x < 0.01 [2, 3, 4]. The saturation hypothesis has been successfully applied to 7 — 7 scattering 
in Ref. [5]. In our recent paper [6] we showed that QCD-based saturation models can lead to 
satisfactory descriptions of cross sections in soft hadron-hadron interactions. 

A pioneering attempt to apply the dipole picture and the physics of high parton densities to 
photoproduction (Q 2 = 0) processes was made by Golec-Biernat and Wusthoff (GBW) within the 
framework of their saturation model [2]. Our goal is to investigate futher to what precision high 
energy photoproduction can be described using the saturation hypothesis. We use a saturation 
model based on a solution of the nonlinear evolution equation [7] . We also extend our analysis to 
the case of nuclei targets. 

The hypothesis of saturation refers to the interactions between partons from different cascades, 
which are omitted in the linear evolution equations (DGLAP and BFKL), and which become more 
significant with increasing energy. The parton saturation phenomenon introduces a characteristic 
momentum scale Q s (x), which is a measure of the density of the saturated gluons. It grows 
rapidly with energy, and it is proportional to ^ [2, 4, 8, 9, 10] with A ~ 0.2. Parton saturation 
effects are expected to be relevant at low values of Q 2 and x, where the parton densities are 
sufficiently large. 

The starting point of this endevour is the successful fit [4] within the framework of QCD, to 
the F 2 structure function data for all values of Q 2 and x < 0.01. This was achieved by using 
an approximate solution to the Balitsky-Kovchegov (BK) [7] nonlinear evolution equation, and 
adding a correcting function to improve the DGLAP behaviour at large Q 2 . Although soft physics 
was not explicitly included, agreement with experiment was found for all quantities associated 
with F2, in particular for the logarithmic slope A = cHni^/c^hi 1/x), where a value of A ~ 0.08 
was obtained for very low x and Q 2 well below 1 GeV 2 , i.e. in the saturation region. This value is 
identical to that of the intercept of the "soft" Pomeron, associated with the Donnachie-Landshoff 
(DL) model [11]. 

In this paper we further persue the hypothesis that colour dipoles are the correct degrees of free- 
dom in photoproduction scattering processes at high energy, even when large transverse distances 
come into play. We adopt the well known expression for the DIS cross-sections 

a^ p (x,Q 2 ) = J d 2 r ± dzP^(Q 2 ,r ± ,z)a dtpole (x,r ± ), (1.1) 

where Q 2 is the virtual photon's four momentum squared. P 7 denotes the probability of finding 
a q — q colour dipole with transverse separation r_|_ inside the photon, s = W 2 is the energy 
squared in the photon-proton system and z, (1 - z) the momentum fraction taken by the quark 
(antiquark) respectively, x = ^ w2 Q + q2) ■ The dipole cross-section a dipole describes the interaction 
of the qq pair with the proton target. The exact form of this cross-section has not yet been 
determined. GBW [2] proposed a very successful parametrization for Odipoie- In our approach 
o 'dipole is obtained by solving the nonlinear evolution equation, which is derived in QCD. 

In [4] the nonlinear evolution equation with certain approximations, was solved numerically, the 
initial conditions were determined by fitting to the low x DIS data. In this paper the resulting 



1 



parametrization for <J dipo i e is adopted for the case of photoproduction. In contrast to our our 
treatment of the DIS case [4], we need to introduce additional nonperturbative parameters to 
insure a finite result in the limit Q 2 — > 0. This leads to a successful description of both the 
photoproduction and low Q 2 DIS data. Our approach for photoproduction is very similar to that 
which we employed for soft hadronic interactions in Ref. [6]. 

As a further test of our hypothesis we apply our model to the case of the photon scattering 
on nuclei targets. This extension does not involve any additional parameters, and can thus 
be regarded as an independent check. All available data on nuclear shadowing is successfully 
reproduced for heavy nuclei. An approach similar to the one presented here for the photon 
scattering on nuclei can be found in Ref. [12]. In [12] a Glauber formula with the GBW model 
was used. In this paper, we rely on the model defined in Ref. [4]. 

The content of the paper is as follows. In Section 2 we discuss the numerical solution of the BK 
equation [7]. We briefly review the GBW saturation model [2], with which we compare some of 
our results. In Section 3 we present the changes that must be introduced to adapt the formalism 
for the calculation of photoproduction cross-sections. Section 4 is devoted to the overall picture, 
including comparison of the model's results with experimental data. We extend our approach to 
nuclei targets in Section 5. Section 6 contains a discussion of our results and our conclusions. 



2 Description of the Saturation Models 

In [4] an approximate solution to the BK non-linear evolution equation [7] was obtained using 
numerical techniques. Below we briefly review the method used and the main results obtained. 
For more details of the method of solution we refer to [4]. 

The solution of the BK equation which we denote by N, takes into account the collective phenom- 
ena of high parton density QCD. Starting from an initial condition which contains free parameters, 
we have numerically solved the nonlinear evolution equation, restricting ourselves to the point 
b± = 0. The parameters which appear in the initial conditions, have been determined by fitting 
to the F 2 data [4], and the resulting approximate solution is displayed in Fig. 1. 
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Figure 1: N(b± = 0) versus r = r±Q s (x). 

The 6j_-dependence of the solution is restored using the ansatz: 

N(r ± ,x;b ± ) = (1 - e -^r ± )s(b ± )^ ^ 
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where k is related to the b± = solution 

K (x,r ± ) = - ln(l - N(r ± ,x,b ± = 0)). (2.2) 
A Gaussian form for the profile function in impact parameter space was taken as 

S(b±) = exp{--^). (2.3) 

71 -Tlproton -^proton 

where Rp roton =3.1 GeV~ 2 refers to the radius of the target proton. Our choice of a Gaussian 
profile function simplifies the numerical calculations. It was shown [13] that calculations per- 
taining to total and integrated cross-sections are insensitive to the details of the profile function, 
unlike the case for differential quantities. 

The dipole-proton cross-section (from eq.(2.4) of [4]) is given by: 

<Jdi P oie(r±,x) = 2 J d 2 b ± N(r ± ,b ± ,x) . (2.4) 
The contribution of the P 2 structure function is 

F 2 (x,Q 2 ) = §^ J d 2 r ± j dz P^{Q 2 ;r ± ,z) cx dipole (r ± , x) . (2.5) 

The physical content of Eq. (2.5) is easily explained. It describes the two stages of DIS [14]. The 
first stage is the decay of a virtual photon into a colorless dipole (qq -pair), the probability of 
this decay is given by P 7 . The second stage is the interaction of the dipole with the target i.e. 
(^dipole in Eq. (2.5)). This equation is a simple manifestation of the fact that color dipoles are 
the correct degrees of freedom in QCD at high energies [15]. The probability P 7 * is given by the 
square of the QED wave functions of the virtual photon, which are well known [15, 16, 17]: 



P^(Q 2 ;r ± ,z) = ^lE^/ x (2-6) 

/ 

{(z 2 + (l-z) 2 )a 2 Kl(ar ± ) + m 2 q K 2 (ar ± ) + AQ 2 z 2 (1 - zf K 2 (ar ± )} , 
with a 2 = Q 2 z(l — z) + m 2 . 

In Ref. [4] we assumed the mass of the quark to be zero i.e. we took m q = in the above 
equation. The structure function P 2 is given by a sum of three contributions: 

P 2 = P 2 + AP + P 2 JV5Q , (2.7) 

where the first term is given by Eq. (2.5). The second term is a DGLAP (large Q 2 ) correction to 
P and is of no importance for the present analysis. These two terms take into account the gluon 
contribution to P 2 . Gluons are related to the singlet quark distributions. The third term in (2.7) 
includes contributions of non-singlet quark distributions: 

F 2 NSQ = E qY . (2.8) 

i=u,d 

At the present stage of our investigation we "borrow" the valence quark distributions (qY) from 
the LO CTEQ6 parametrization [18]. It is important to note that these distributions decrease 
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with decreasing x, and are practically of no significance below x ~ 10~ 3 . As the valence quark 
parametrization is not given below some scale of order lGeV 2 , for low photon virtualities we 
freeze the valence quark distribution at its Q 2 = 1 GeV 2 value. 

A popular saturation model was proposed by GBW [2]. In this model the dipole cross-section is 
taken as 



&dipole(r_L,x) — (T 



r 2 . 



1 - e ^(-77^) 



(2.9) 



ar 2 

, with Rl(x)[GeV~ 2 ] = (^) A - The values of the parameters which were determined by fitting 
to DIS data at HERA for x < 0.01, are: a = 23 mb, A = 0.29 and i = 3x 10~ 4 . The 
rj_ dependence was taken as Gaussian, which leads to a constant-cross section a for large r\ (or 
small Q 2 ). 

The dipole cross-section of the GBW saturation model and N are closely related. Both models 
take into account gluon saturation, preserve unitarity, and describe the physics associated with 
"long distances". Whereas N has its origins in QCD, (Jdipoie is phenomenological. 

There are differences between the two models. Unlike the GBW <7di po ie, the dipole cross-section 
obtained from the solution of the BK equation is not saturated as a function of x. With the 
assumed Gaussian profile function, the integration over impact parameter b± leads to a logarithmic 
growth with decreasing x. 



3 Details of Our Calculation 

In Ref. [4] we considered a massless (m q = 0) case only. To proceed to the limit Q 2 — > we need 
to introduce a finite mass as a cutoff for the r± integration in Eq. (2.5). It is important to stress, 
however, that we do not modify the dipole cross section <Jdi po ie (2.4), which was constructed in [4] 
from a fit to the Fi data. Since x is not defined for photoproduction, to relate x to the energy 
of the process we need to introduce an additional nonperturbative scale, denoted as Q\. Though 
these two parameters are not obviously related, we set Ql = Am 2 . We redefine x as 

x = (Q 2 + Ql)/W 2 ■ s = W 2 , (3.1) 

and use m q as a fitting parameter for high energy photoproduction. The energy dependence of the 
photoproduction cross section enters only through the ^-dependence of the dipole cross section, 
the later being adjusted or constructed to describe DIS data of the F 2 structure function. 

Note that having introduced these modifications we alter the results in low Q 2 domain. This can 
potentially spoil our original fit [4]. Thus photoproduction should be considered simultaneously 
with DIS. The fit to the later should not be affected by the introduction of new parameters. 

To obtain a finite result in the limit Q 2 — > 0, GBW [2] introduce the quark mass as a regulator, 
and assume a common mass of 140 MeV for the three light quarks, "which leads to a reasonable 
prediction in the photoproduction region" . We refer the reader to [2] for details of the GBW model, 
in this paper we use the GBW model only to compare with our results for photoproduction on 
nuclear targets. 

In the colour dipole picture, one can hopefully reproduce only the asymptotic energy dependence 
i.e. the Pomeron contribution. The solution to the BK equation is also used to model this 
vacuum exchange channel in the nonperturbative domain. Thus parton-hadron duality is used in 
this approach. At lower energies where most of the photoproduction data are available, there are 
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also contributions from nonvacuum exchanges. The evaluation of these contributions is beyond 
the scope of our model. 

In [4] the contribution of the nonsinglet exchange was represented by the valence quarks. Since 
the latter were frozen at Q 2 = 1 GeV 2 , their contribution to the cross section would have the 
wrong dependence at Q 2 — > (become infinite). An improved form for this contribution of the 
valence quarks below 1 GeV 2 is: 

f nsq _ (Q 2 ) 1+l3 (IGeV 2 + /i 2 ) v 2 y 2 _ 2 

t2 - 1 (GeV 2 )^ Q 2 + p 2 ^ d 6l qAQ ~ 1 GeV } ■ {6 - 2) 

Vector current conservation requires (3 > 0, in principle, both f3 and \x are fitting parameters. 
This freedom was not used in [4] where f3 — and fj 2 — > 0. The energy dependence of the valence 
quark contribution enters through its x-dependence, which is frozen at Q 2 = 1 (GeV 2 ). Taking 
fj, as a fitting parameter produces a good description of both the DIS and photoproduction data 
(keeping (3 = 0). 

In our approach we use parton-hadron duality in two ways. For the vacuum exchange we assume 
that the saturation result, which was derived for large Q 2 virtualities, is a valid description (on 
average) of the cross section at small values of Q 2 . For the contribution of the non-singlet structure 
function we can use a parametrization that is successful for soft interactions, namely the exchange 
of secondary Reggeons. We believe that this exchange can be used for larger Q 2 (up to 1 GeV 2 ). 
In the spirit of Ref. [11, 19] and using Gribov's formula [14] the secondary Reggeon contribution 
is 

s\ Q « rdM 2 p(M 2 ) . o \ , n „v 

7j /wm^ r(M) ' (3 ' 3) 

where p(M 2 ) is related to the cross section for e + e~ annihilation while gn is a residue of the 
secondary trajectory. We approximate (3.3) by the following contribution 

M 2 

a" ~ /(0) — - (s/s o r« ; a R = -0.45 ; s = 1 (GeV 2 ) . (3.4) 

f(0) is the residue at Q 2 = 0, and is determined by fitting to the low energy photoproduction 
data. M is an additional fitted parameter. 



4 Photoproduction and DIS at low Q 2 

A fit to high energy photoproduction data yields the value of 

m q = 0.15 GeV. 
Note that we obtain practically the same quark mass as GBW [2]. 

As most of the photoproduction data was measured at lower energies, it is essential to include 
an additional component which is dominant in this energy range. We attempt two alternative 
parametrizations for this component. 

The first is based on the Regge pole model Eq.(3.4) which we denote: 

• Model ST: In this model there are two parameters to be fitted which determine the 
contribution of the non-leading (or secondary trajectories) to the DIS cross-section (see 
Eq. (3.4)). These are 

/(0) = 0. 19 mbarn; M 2 = 2 GeV 2 . 
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The cross section obtained with these parameter values is indicated by the dashed curve in 
Figs. 2, 3. The second parametrization is based on the quark model: 

• Model VQ: This alternate method of accounting for the non-singlet contribution to the DIS 
cross-section is via the valence quarks (see Eq. (3.2)). There are also two free parameters 
here, whose values are determined by fitting to the data. These are 

(3 = (kept fixed); fi 2 = 0.13 GeV 2 . 

The resulting photoproduction cross section with this parametrization for the valence quarks 
is shown as the full line in Figs. 2, 3. 



200 
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Figure 2: Photoproduction as a function of 7 — p energy. Preliminary data from [20]. The solid 
line indicates the results of Model VQ, while the dashed line those of Model ST. 

We note an interesting duality between the two models, both produce reasonable descriptions 
of the experimental measurements (see Fig. 2), including the low energy photoproduction data. 
For the case of a proton target we do not present a comparison with the GBW model, since the 
later was formulated only for high energies and does not include any contribution of secondary 
trajectories. 
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Figure 3: DIS cross section at very low Q 2 The two curves correspond to two different models for 
the secondary trajectory (see text). Data from [21]. 



7 



5 Photoproduction on Nuclei 



In our approach we obtain a very small effective radius of the nucleon. The consequence of such 
a small radius is that the typical parameter that governs the dipole interaction with the nuclear 
target, is also very small, namely 



71 ^proton 5^(0) 



< 1 



for A < 200, 



(5.1; 



where SA_{b±) denotes the Wood-Saxon impact parameter distribution of the nucleons in a nucleus. 
In the standard approach, 7ri?L. oton Sa(0) is considered to be large. Instead of one equation [7], 
Eq. (5.1) leads to a set of non -linear equations which have to be formulated and solved. Before 
doing so it is necessary to determine the b± distribution for dipole-nucleon scattering, which is 
still an ansatz in our approach. At present, the best that we can do is to extend the Glauber 
formalism which takes into account all rescatterings inside the nucleus. This approach is certainly 
correct at not too high energies, however to determine the region of applicability we need to solve 
the explicit set of equations. 

The Glauber formula has a familiar form: 



a(-f*A — >j*A) = / d 2 b ± / dzd 2 r ± P^(r ± ,z,Q 2 )2 



(5.2) 



In our model (Jdipoie is given by Eq. (2.4). We use the GBW model for the dipole cross section 
for comparison (see also Ref. [12] for a similar analysis). 

As discussed previously, at low energies we need to add contributions of the secondary trajectories. 
For the Glauber approach the A dependence of the secondary trajectory can be estimated as 
follows: 



d 2 r x \^(r x )\ 2 d 2 b ± S A (b ± )e 



y dipole( r ±' x '> n \ 
— o bA(0±) 



(5.3) 



Here denotes the pion wave function. For ty n we use a simple Gaussian parametrization taken 
from Ref. [22]. The physical meaning and the derivation of Eq. (5.3) is clear from Fig. 4. Indeed, 
the exchange of secondary Regge trajectories is screened by the black disc type interaction in the 
vacuum channel (the two diagrams in Fig. 4 which cancel each other at small b±). The factor 
exp (— <JdipoieS(bj_) /2) is the probability that the dipole does not to have any inelastic interaction 
(for a more detailed discussion see Ref. [23]). 




A 



- r: 



N 



,GM 



Figure 4: The diagrams for the secondary reggeon interaction with nuclei. N^ M is the dipole- 
nucleus amplitude in the Glauber-Mueller approach. The two diagrams cancel each other and 
lead to Eq. (5.3). 

We compute the total cross section for 7 — A collision as a sum of two terms 



a^*A) = Eq. (5.2) + Af ff x Eq. (3.4). 



(5.4) 
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Figs. 5, 6 7 illustrate the dependence of A e ff = a{^*A)/a{^*p) as a function of x for different 
values of the photon virtuality Q 2 for photo scattering on Pb, Xe, and Ca respectively. Fig. 
8 shows a comparison for F 2 on Ca. The data sets have a few percents overall normalization 
uncertainity. 

In the figures below we use the following notation: the solid line shows the results of our formalism, 
where <Jdi P oie m Eq. (5.2) is calculated from the solution of the non-linear equation. The dotted line 
illustrates the Glauber approximation for the GBW model for the dipole-nucleon cross section, 
without a secondary Regge pole contribution. 1 

Agreement with experimental data for both dipole models is good. A comparison of our model 
prediction for a lead target with the results other models can be found in the recent preprint [24]. 

In Fig. 9 we plot the x dependence of a for the parametrization a(j* A) oc A a . Compared to 
GBW our results appear to be in slightly better agreement with the data. However, the data are 
not precise enough to really distinguish between the two models. 

Our prediction for the Q 2 dependence of a is shown in Fig. 10. 

We compare the predictions for the energy behaviour of our approach with the relevant data for 
A e ff/A for photoproduction on Pb in Fig. 11, and on Cu in Fig. 12. We slightly underestimate 
the data for Pb case, and are consistent with the Cu measurements. Some points of the Caldwell 
data (1979) [28] are read off the plots and divided by our photoproduction cross section. 

The curves in Fig. 13 which show predictions for the values of the total 7*A cross section are 
crucial as a check of our approach. The results of both models are very close and this boosts our 
confidence in our calculations. The predictions for F 2 for 7* - gold scattering for both approaches 
are given in Fig. 14, again the results are close to one another. 

6 Discussion and Conclusions 

This paper demonstrates that assuming QCD saturation of the parton densities at high energies, 
one is able to describe photoproduction processes, as well as DIS at small values of photon 
virtuality. We establish that the revised parametrization does not spoil the predictions made for 
DIS processes at large and moderate values of Q 2 . This is important as it allows us to use the 
DIS data to fix the phenomenological parameters of our approach. 

As in our treatment of the soft hadron interactions in Ref. [6], the present analysis is based on 
two parameters and two major assumptions. The parameters are: 

(a) The effective radius of proton which is taken from the description of the DIS data for x < 1CT 2 ; 

(b) The parameter Q%, (introduced in Eq. (3.1)), is adjusted by comparing with the high energy 
photoproduction data. 

The two assumptions are: 

(i) We use Eq. (2.1) as an ansatz for the 6j_-dependence. This Glauber- type b± dependence has 
no solid theoretical foundation and should be considered to be an initial guess. It is therefore 
reassuring to see that our successful description of the photoproduction data is in reasonable 
agreement with this ansatz; 

(ii) We use hadron-parton duality to introduce the contribution of the secondary Regge trajecto- 
ries. The result that the contribution of the secondary reggeons is dual to the valence quark one, 
is unexpected and interesting. 

1 In principle, contributions of the secondary trajectories should be added to the GBW model as well. 
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Figure 5: A e ff/A on a Pb target. The solid line denotes the results of our model including a 
secondary Regge contribution. The dotted line is the result for the Glauber approx. for GBW 
model without a secondary Regge contribution. Data from [25]. 
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Figure 6: A e ff/A on a Xe target. Notation as in Fig. 5. Data from [26]. 
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Figure 7: A e ff/A on a Ca target. Data from [27]. 
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Figure 8: F 2 on a Ca target. Data from [27]. 
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Figure 9: a as a function of x for different values of Q 2 . Notation as in Fig. 5. Data from [25]. 



14 



0.9^ 


i 


III 1 1 III 


mi i i i i i i i i 


0.92 






/ / 


3 0.9 


x = 
x = 


10-3 ^^-^^ 

10-4 ^ — 




0.88 


x = 


10-5 




0.86 


1 


i i iii 


iii i 



lO- 1 10° 10 1 



Q 2 (GeV 2 ) 

Figure 10: a as a function of Q 2 for fixed values of x. Prediction of our model including a 
secondary Regge contribution. 
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Figure 11: Photoproduction: A e ff/A as a function of 7 energy for scattering on a Pb target. 
Data from [28]. 
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Figure 12: Photoproduction: A e ff/A as a function of 7 energy scattering on Cu target. Data 
from [281. 
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Figure 13: The predictions of our model for the DIS cross section on Au at very low Q 2 . Notation 
as in Fig. 5. 
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Figure 14: F 2 for gold (A = 200) as a function of x for fixed values of Q 2 . Notation as in Fig. 5 
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Our main result, that we are successful in describing the photoproduction data for nuclear targets, 
is compelling and encourages us to suggest further measurements of photon - nuclei interactions. 
We stress that the Glauber approach for nuclear targets does not contain any additional assump- 
tions, thus nuclear data can be used for clarifying the situation regarding our assumptions and 
fitting parameters. 

This paper and our previous work [6], provide evidence that an approach based on the QCD dipole 
picture and on QCD saturation is successful in describing processes which traditionally have been 
associated with soft physics. They also reinforce the results obtained in Ref. [4] namely, that 
QCD saturation is a natural bridge between soft physics and conventional perturbative QCD. 
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